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Abstract 



The HarishChandra-Itzykson-Zuber integral over the unitary group U(k) 
((3 = 2) is present in numerous problems involving Hermitian random ma- 
trices. It is well known that the result is semi-classically exact. This simple 
result does not extend to other symmetry groups, such as the symplectic or 
orthogonal groups. In this article the analysis of this integral is extended 
first to the symplectic group Sp(k) (/3=4). There the semi-classical approx- 
imation has to be corrected by a WKB expansion. It turns out that this 
expansion stops after a finite number of terms ; in other words the WKB ap- 
proximation is corrected by a polynomial in the appropriate variables. The 
analysis is based upon new solutions to the heat kernel differential equation. 
We have also investigated arbitrary values of the parameter (3, which char- 
acterizes the symmetry group. Closed formulae are derived for arbitrary (3 
and k — 3, and also for large (3 and arbitrary k. 
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1 Introduction 



The HarishChandra and Itzykson-Zuber (HIZ) integral |], ^ |Jfor the unitary 
group Uk is defined by 

I = J dtfexpitrU+XUA (1) 

where X and A are diagonal k x k Hermitian matrices. It is well known that 
this integral is given, up to some normalization of the measure, by 

A(X)A(A) 1 j 

where A(X) and A (A) are the Vandermonde determinants of the eigenvalues 
of X and A, (A(X) = J^(xi — Xj)). This HIZ formula may be easily derived 
by considering the Laplacian operator pj H] 



d 2 

Its eigenfunctions are plane waves 

Le itrAX = ( trA 2) e itrAX ^ 

One can construct a unitary invariant eigenfunction of L, for the same eigen- 
function trA 2 , by the superposition 

df/e 4trAUXU \ (5) 



which is nothing but the HIZ integral. The integral / beeing unitary invari- 
ant, is a function of the k eigenvalues Xi of X. The one dimensional quantum 
mechanics of the eigenvalues, for unitary invariant eigenstates, is then easily 
solved in terms of free fermions, leading to the unitary HIZ formula. 

The same considerations hold for the three ensembles (3 — 1, 2 and 4, 
corresponding to the orthogonal, unitary and symplectic ensembles, with 



The Laplacian acting on invariant states, may be expressed in terms of a 
differential operator on the eigenvalues : 

k f)2 k id 

+ £ ^^^ / = - eJ ' ( 7 ) 

i=l i=1) (j^) a j A j 



with the eigenvalue e 



e = 

8=1 



Note that the integral / is manifestly symmetric under interchange of the 
matrices A and X, but the procedure is dissymetric. The solutions will of 
course restore this property, which is not obvious if one considers the equation 
(0) alone. 

The x-dependent eigenfunctions of this Schrodinger operator have a scalar 
product given by the measure 

{VilVz} = dxi---dx k \A(xi---x k )\P ip* 1 (x 1 ---x k )<P2(xi---Xk) (9) 



The measure becomes trivial if one multiplies the wave function by lAj' 3 / 2 . 
Thus if, for some given ordering of the x^s, one changes I(x) to 

^(xi-.-a*) = \A(x 1 ---x k )f /2 I(x 1 ---x k ), (10) 

one obtains the Hamiltonian, 

E |? - - !) E T^hr^ = -«/>■ (") 

i=l ^ i<j l x « 

The relation between matrix quantum mechanics and many-body problems 
with 1/r 2 pair potentials was already present in |7j and it has also be used 
for the study of Selberg integrals by Forrester |8j . This Schrodinger equation 
is a simple Calogero-Moser model. 

Although the solutions of this many-body problem are known, they are of 
little use for our purpose, since it is already quite involved to recover from the 
formal expression of the solution for arbitrary (3, the simple unitary result of 



HIZ and it looks formidable to extend it to other values of /3 []T0|]. For (3 = 2, 
the solution is again given by plane waves in the Xj and (taking into account 
the symmetry under permutations of I), one obtains the HIZ formula. For 
(3 = 4, the problem is much less trivial, but it leads to explicit solutions for 
finite values of k. Indeed the problem turns out to possess simple solutions 
of the form of symmetrized sums of plane waves multiplied by polynomials 
in the variables — Xj)(Xi — Xj), providing therefore a complete explicit 
solution for the Sp(k) group integration. 

For general k, the solution of (|TT|) is of the form 

^ = e *(Ai*i+-+A** fc ) x ( 12 ) 

where x satisfies 

E ^ + 2 * E ^ - y E z^p]* = ( 13 ) 



in which 



y = /?(§-!). (14) 



For k=2, it is immediate to verify that there is a solution, function of the 
signle variable r = (Ai — \2}( x i ~ x %), whose asymptotic expansion for large 
t is 

x = i+EH/rrn(p-i-|-) 

n=l 1 Z ^ 



= I u (z)V2nze~ z (15) 

in which the Bessel function I v {z) is defined with v 2 = (y + l/2)/4 and 
z = ir/4. 

When /? = 4, the operator ^ E 7 V2 annm ilates the function 

j=l j<j \ X i ~ x j) 

A _1 (xi, • • • , Xk) ■ Consequently the solution of ( fL3l) may be written 

X(xi • • • x fe ) = — r (16) 

A(xi • • - xjfc) 



in which f{x\ ■ ■ -Xk) is a polynomial of degree k(k — l)/2 in the x a 's. Defining 

Tij = (Ai- \j){xi-Xj) (17) 

one finds for k = 3 



, 1 1 1 , , 1 1 1 s 
X = 1 + 2i{— + — + —)- 4( + + ) 

712 T 23 T 3 l r 12 T 2 3 T23T31 T31T12 

-12i . (18) 

T-12T23T31 

Remaining still with (3 — 4, one may factor out the Vandermonde deter- 



minant from x as (16), and one obtains a differential equation for /, 



g + 4 - 1) + *V)(A--) + 2i g ^ = (19) 

The solution of this equation is obtained by expanding in powers of the Aj's, 
but the expansion ends with the highest degree, the Vandermonde product 
of degree six A(Ai • • • A 4 ). Using the notation of Q17D , = (xj — xj)(Xi — 
we obtain 

/ = 1 - l -{ri2 + Tig + T U + r 23 + r 24 + T34) + 0(r 2 ) (20) 

as shown in the appendix of 0. In the next section, we consider the sym- 
plectic case {(3 = 4) for larger values of k, and we will clarify the meaning of 
the coefficients which appear in the expansion of x- 

2 Decomposition into complete graphs in the 
symplectic case 

Let us discuss now the asymptotic expansion of x f° r P — 4, i.e. for the 
symplectic group integration of the HIZ integral. Such integrals appear in 
various problems ; for instance one encounters it in the calculation of the 
average of products of characteristic polynomials of random matrices in the 



Gaussian orthogonal ensemble < Ili=i det(Aj — X) >. Indeed we have shown 

in a previous work ||] that this average may be expressed as an integral over 

a k x k quaternion matrix ; this representation exhibits a duality between the 

GOE average and a symplectic HIZ integral. Thus this requires to compute 

the HIZ integral for (3 = 4, as discussed in , in which the expressions 

were given up to k = 4. We extend here the construction of the polynomial 

solutions to arbitrary values of k . 

As discussed before, the expansion of the function x m inverse powers of 

the variables l/r^, defined in fllTD , terminates with a term proportional to 
1 

— , of degree k(k — l)/2. The monomials in the variables l/r^ may be 

i<j Ti 3 

represented graphically : one marks k points, and a line between the points 
i and j will represent a factor — in a given monomial. A complete graph is a 

Tij 

graph in which all the pairs of points are connected by a line. For instance, 
a triangle (k=3) and a tetrahedron (k=4) are complete graphs. The k-point 
complete graph has k(k-l)/2 lines. A complete graph of k points receives a 
factor Ck in which 

k 

C k = f[l\ (21) 

i=i 

This number is related to the integral, 

= J^JU^-z^e-^lldz, (22) 
The last term of highest degree in x , is equal to 

' k C k \{— (23) 



i 



i<j r n 



All the terms of the expansion of % have coefficients which may be decom- 
posed into combinations of factors C n with n < k. 
For instance , the k=3 solution is 

C C (C ) 2 

X = 1 + i—^r-{v\ +v 2 + v 3 ) + i 2 — ^—(viv 2 + v 2 v 3 + v 1 v 3 ) +i 3 C 3 viv 2 v 3 (24) 
Go Oi 

5 



where 1*3 = — , V\ = — ,and v 2 = — . The graphs of the various terms in \ 

Tl2 T 2 3 T 31 

may all be analyzed as a decomposition into complete graphs. For instance, 
V1V2 consists of two lines connecting the point number three to points one 
and two. The lines (2-3) and (3-1) are complete graphs of two points, hence 
the factor Cf, with the point three counted twice, thus a division by (Ci), 
equal to one. (One has to multiply further the result by a power of i raised 
to the number of lines) . 

When k = 4, by this rule of decomposition of the graphs into complete 
graphs, we find the coefficients of all the terms in the form 

ri ri 6 
-. . .^2^1 

°o i=1 

(C ) 2 

+i 2 2 (V1V2 + V 2 V 3 + VtVs + + ^1^5 + V2V5 + V 2 Vq + V 3 Vq 

+v 3 v A + v 4 v 5 + v 4 v 6 + v 5 v 6 ) 
+i 2 (C 2 ) 2 (v 1 v 6 + v 2 v 4 + v 3 v 5 ) 

~77^ (viV 2 v 3 + + v 2 v 5 v e + v 3 v 4 v e ) 

(C 2 ) 3 

77^(^1 ^3^5 + ViV 3 Vq + V 2 V 3 V 4 + V3V4V5 + V 3 V b V & + ViV^Vq + V 2 V 4 V<s 
+V 2 V 3 V<j + ViV 2 V A + ViV 4 Vq + ViV 2 V 6 + V 2 V 4 V 6 ) 
+i 3 2 (ViV 2 V 5 + ViV 3 V 4 + V 2 V 3 V 6 + V 4 V 5 V e ) 



(C 2 ) 4 C , 
(C1) 4 



+i in \A ( V 1 V 2V4V 6 + V!V 3 V 5 Vq + V 2 V 3 V 4 V 5 ) 



C 3 C 2 

n {v 2 V 4 V<zV§ + V 2 V 3 V h V§ + VxV^Vz + V 3 V 4 V^ + V X V 3 V 4 V§ + V 2 V 3 V 4 V§ 

w 

+ViV 2 V 3 V b + ViV 2 V 3 V§ + ViV 2 V 3 V A + ViV 3 V A V b + ViV 2 V 4 V b + ViV 4 V b V & ) 

(C3) 2 

^2 

+V 1 V 2 V 3 V 4 V 6 + ViV 2 V 3 V 4 V 5 ) 

+i <i C 4 v 1 v 2 v 3 v 4 v b v & (25) 
where we have used v 1 = l/r 12 ,v 2 = l/r 23 ,v 3 = l/r 13 ,v 4 = l/r l4 ,v 5 = 



6 



l/ r 24, vq = l/r 34 . Again the graph of the last term V1V2V3V4V5V6 is a complete 
graph of 4-points. Substituting in this formula the values C n = H"=i j!, one 
recovers precisely the k = 4 result, given in an earlier publication [f|. ||. 

The decomposition of graphs into complete components has certain gen- 
eral structures for arbitrary values of k. The term of highest degree is 
jk{k-x)/2(j k . nex ^ on g ; obtained by deleting one line, has a decompo- 
sition given by 

M (26) 

For k=5, this decomposition means that the graph is made of two tetrahedra 

sharing a triangular face ; hence the division by C3 to avoid double counting 

of this face. This explains the above factor fl26|). In the case k=4, this factor 

(C3) 2 

represents two triangles connected by one edge, and it reads — — — = 72. 

C2 

When two lines are eliminated from the maximum complete graph C^, 
we have two type graphs. 

(i) The first case is the elimination of two lines emerging from the same 
point. In this case the decomposition is Ck-iCk~2/Ck~3- For k=4, it is a 
triangle and one line connected at the vertex. For k=5, it is a tetrahedron 
connected to a triangle through one edge. 

(ii) The second type consists of deleting two non consecutive lines. In 
this case, we have the decomposition of (Cfc_ 2 ) 4 Cfc_4/ (Ck-3) 4 ■ In the case 
k=4, it corresponds just to a square, four line connecting the four vertices of 
a square. For k=5, it is a graph of four triangles connected to each other by 
four edges. 

If we divide this factor (Cfc_ 2 ) Cfc_4/ '(Ck-3) 4 by C/., as appeared in the 

13 2 

expression for / in ft becomes — , — , — for k= 4,5 and 6, respectively. 

18 80 75 

These values have been found by a direct evaluation of the solution of the 
differential equation ([19]). Thus we have, for (3 = 4 and for arbitrary values 
of k, an expression for x, whose coefficients are given by the decomposition 



7 



of the complete graphs C n . 

The case (3 = 4 is thus easy to analyze for two reasons. First there are no 
"multiple" lines in a graph, meaning that every monomial in \ is °f degree 
at most one in any variable 1/%. Then the factor C n = Ylj\ for complete 
graphs has a simple structure. However, for (3 ^ 4, the number of lines 
between two points may be arbitrarily large. Then the contribution C n of 
complete graphs has a more complicated expression. For values of (3 of the 
form 

f3 = 2(m + 1) (27) 

a complete graph between two points consists of m lines joining these two 
points and 

C 2 = h-l + h(-2 + -)■•■ (-(m - 1) + — ) (28) 
2 V 4 A 6 J V V ' 1m' y ' 

in which y has be defined in fll4l) . The evaluation of C3 is non trivial and 
we will discuss it in the following section. For /3's of the form ( p7[ ) the 
previous graphical analysis of decomposition into complete subgraphs, may 
be extended to arbitrary k. However, even then, more complex rules for non- 
complete graphs are required. To understand this complexity, we consider 
the k = 3 case for arbitrary (3 in the next section. 

3 k=3 for arbitrary (3 

Even for the simple k = 3 case, the HIZ integral is not known explicitly 
in the form of a WKB expansion. An analysis in terms of the solutions of 



the Calogero three-body problem has been performed in fi~0 |. However, its 
outcome is far from our purpose since it is already quite involved to recover 
from there the unitary HIZ formula, and one does not see that for values of (3 
such as four, the WKB expansion is a polynomial. Therefore we start again 
with the partial derivative equation ([13]) for k = 3, and expand the solution 



8 



as 



y z n + m+r C n , m , r 
n,m,r=0 r 12 r 23 T 31 



Note that it is not at all obvious that the equation (|T^) admits such a solution. 
Indeed, assuming that the solution is only a function of the r's, one finds 
easily that 

I* E A ^-E TZTZ^U = "(Ax - A 2 ) 2 r? 2 [2i^- + y] x + perm. (30) 

i=l 2 i<j V 1 j) 12 

However the term involving the second derivative does not share this struc- 
ture since 

3 «2„ Q2„ 



2(Al " A2) [Tl2 ^ + Tl2 ^ J 



-2(Ai - A 2 )(A 3 - AxKVi + perm (3i) 

The last term, with cross-products in the A a 's, is a priori different. However 
there are, order by order, a set of identities which allow one to cast those 
cross-products into sums of squares. Those identities follow from {x\ — x^) + 
(x 2 — £3) + (X3 — xx) = 0, leading to 

(Ai - A 2 )r 2 3T 3 i + perm = 0. (32) 

From there one derives order by order the identities which are needed to cast 
the cross-products into squares. Let us quote the simplest 

(Ai - A 2 )(A 3 - Ai)r 2 2 3 + perm = (X x - A 2 ) 2 r 3 ir 23 + perm. (33) 

Order by order those identities allow us to satisfy the conditions on the 
expansion of x, which make it a function of the r variables . One finds easily 
that C n oo, which is associated to the diagram consisting of n lines between 
two points, is given by : 



n 



C n00 = II(-(i - 1) + l) (34) 
1=1 u 



The generating function of those coefficients is the modified Bessel function, 
solution of the k=2 problem. Note that C„ jmjr is symmetric under permuta- 
tions of n,m and r ; for instance C mj0 ,o = Co,m,o- 

When r=0, one finds that C n ^ m $ has a decomposition into a product of 
two factors, 

Cn,m,0 — Cn,0,oCm,0,0- (35) 

For the full HIZ integral, triangle graphs, in which the three integers 
(n, m, r) of C n , m ^ r are non-zero, play an essential role. 

From the differential equation one can derive a recursion equation 
for the C n , m ^ r . However a direct expansion of the solution of is extremely 
tedious and complicated at increasing orders, since one has to use more and 
more identities of the type (|53"D, but of higher degree. One needs triagle 
identities. The expansion for the differential equation has a simpler 
structure ; it gives a relation between C n>mjn with n + m + r = P + 1, in 
terms of C n ^ m ^ r with n + m + r = P ; this relation is valid for a fixed given 
value of m. If one writes C n , m ^ r as a linear combination of C n >^ m y, where 
r' < r, and n' + r' = n + r — 1, one obtains 



C„, m ,i = C n , mfl {nm + |) (36) 

Cn,m,2 — Cn,m,l( — — 1 + ~) + C n+ i tln fi — (n + 1) (37) 

C nm 3 = C nm 2{ — { n ~ 2) — 2 + — ) + C n+ i m i— — 

3 o 3 

+C n+2tmfl -{n + 2) (38) 

For general r (r > 3), we have 

m(n-r + l) ^_ 
r 2r 

r j I I ] / 

+C n ^ r -2 — {n - r + 3) + Cn^m.rs — (u - T + 5) 
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+ ••• 

Tfl 771 

+C n+r _ 2 , m ,i — (n + r - 3) + C n+r _i >m;0 — (n + r - 1). (39) 

The solution of these recursion equations guarantees automatically the sym- 
metry under exchange between n, m and r : C n>m>r = C n ^ m = C m ^ r = 
n _ n - n 

^m,r,n ^r^jTU ^ J r,m,n' 

From these recursion formulae, one may determine iteratively all the co- 
efficients Cn^m^. The resulting expressions are complicated for general y. For 
instance, even at low order, one finds 

^ 2 , 2 =(f) 2 (-l + |) 2 (-6 + ^+^) (40) 

^3,3,3 = (|) 2 (-1 + V -n-2 + |) 2 (-51 -fy + + ^) (41) 
The solution for (3 — 4 and A; = 3, may be written in the compact form 

X = ^ (42) 
= 2i(— + — + — ) - At (43) 

r 12 r 23 r 13 r 12 7 "23 r 13 

in which it is understood that one expands it in powers of 1/r, dropping 
any term involving powers of 1/Vjj greater than one. This corresponds to the 
absence of multiple bonds between two points in the graphical representation 
for (3 = 4. This result agrees with (|18D. This certainly suggests that ought 
to be written in terms of Grassmann variables. In the next section, we 
consider such expressions for general y in an expansion in 1/y. This also 
provides consistency checks for the recursion equation (|39f) . 

4 Large y expansion 

One can obtain an expansion in powers of 1/y if one returns to the differential 
equation (|13D, taking now as variables 

u a = x a /y. (44) 
11 



The function x( u i ' ' ' u k) is solution of 



l 1 -£^+*£^i r -T,7^-^}x = 0- (45) 
The leading term xo is thus solution of 

[2i S<-£CT ] ^ a (46) 

We look for a solution which is a symmetric function of the variables 

i 

v a,b = tt t^t r , a < b (47) 

(A a - M){u a - u h ) 



which thus satisfies 



The equation ( f46| ) takes then the simple form 

E (A a -A 6 ) 2 ^[2-^-l] Xo = 0. (49) 

\<a<b<k OVab 

There is a trivial symmetric solution which satisfies the k(k — l)/2 equations 

[2^- - l]xo = 0, (50) 



namely 



Xo = exp (i X>) = exp (f g (A n — X b )(x a — x b ) ^' 



One thus writes 



X = ^ 

1 1 

4> = o + -0i + ^02 + ---, (52) 

y y 2 

with 

= \ E ( 53 ) 

Z a<b 
12 



The next terms are solutions of 

? ( W + iw )+2l ? A °*: =0 ' (54) 
? (2 ^^ + ^ )+2! p^ = a (55) 

In order to make it clear that it is a non trivial property of these equations 
to have solutions which are symmetric polynomials in the variables v a b, we 
shall examine these equations in some detail. 

For simplicity of notations we shall limit ourselves to k=3, the general- 
ization to arbitrary k beeing immediate. It is thus convenient to simplify 
slightly the notations and write 

viz = V-i , v 23 = V\ , U31 = V 2 (56) 



Let us first consider ( p4| ) ; since the function <p\ is a symmetric function of 
the Vi's one finds easily, as in that 



E = i[(A 2 - A 3 ) + perm.]^ (57) 



The explicit solution (p3[) allows one to write 

E((|^) 2 = -^[(A 2 -A 3 )V 1 4 + perm.] + i[(A 1 -A 2 )(A 3 -A 1 )V^ + perm. 

= -2[(A 3 -AsW + penn.] (5* 

The second term of the first line of (|57j) does not involve a priori (A a — Ah) 2 
; however the variables Vi are not independent since they are defined in terms 
of the differences (ti — t 2 ), (t 2 — £ 3 ), (t 3 — ti) whose sum is zero. Thus 

V 2 V 3 {\ 1 - A 2 )(Ai - A 3 ) + perm. = (59) 

and one can check without difficulty that 

VfV?(Ai - A 2 )(Ai - A 3 ) + perm. = -ViV 2 V 3 [(A 2 - A 3 ) 2 Vi + perm.]. (60) 

13 



With the help of this identity the equation for 0i takes the simple form 

(A 2 - As) V x 2 [^ + -Vl + V x - ]v 2 V 3 ] + perm. = 0. (61) 
OV\ 4 4 

Thanks to the identity floTf ) one can look for a solution of the three equa- 
tions 

T + \ V i + V i~ Va = (62) 



914 4 1 1 4 
plus the other two deduced from the cyclic permutations. The system is 
trivially integrable and one finds 

h = \v x v 2 Vz - hyl + v 2 2 + v 3 2 ) - ^(n 3 + v 2 3 + vi). (63) 

For arbitrary k this generalizes to 

01 = 1 Y VabVbcVca ~ £ Y V ab ~ ^ Y V lb, (64) 
(abc) a<b a<b 

in which the sum (abc) runs over all distinct triplets of indices. 

We have exposed this simple calculation in some detail, in order to make 
it clear that the success of the method relies on two successive facts : 



• Identities such as |6"UD which allows one to write the differential 
equation as three separate equations for d(f>/dV a 

• The integrability conditions of this system of three equations are sat- 
isfied. 

At higher orders the technique is identical. It requires new identities of higher 
degree, such as 

V?Vi(Vi + If) (Ai - A 2 )(A! - A 3 ) + perm. 

= -(A 2 - A 3 ) Vi^^ift 2 + V 2 2 + V?) + ^V 1 V 2 V 3 \ + perm.. (65) 
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and many others ; thereby one obtains again remarkably an integrable system 
of three equations. At order 1/y 2 the solution is 

02 = (^i 3 + vi + vi) + \{v? + v 2 A + vi) - \v 1 v 2 v^{v l + v 2 + v 3 ) 

+^(V? + V 2 5 + V 3 5 ) - + V 2 2 + V 3 2 ) (66) 

5 A cubic identity 

Let us consider the case of k = 4 for different values of /3, assuming that it 
is an even integer. For the tetrahedron (k=4), the coefficients of the expa 
nsion of x may be written as 

3 fi 
y _ jn+m+r+p+q+l ^nmr;pql /g^x 

n,m,r,p,q,l=0 r 12 r 13 r 23 r 24 r 34 r 14 



where j = (3/2 — 1. By imposing to satisfy the differential equation (|T9|) , one 
should determine those coefficients. For instance, in the case (3 = 6 (or y = 
12), one finds C 222;222 = 8-7(6-5) 2 (4-3) 3 , C 222;221 = (2-7)(6-5) 2 (4-3) 3 , C 222;21 i = 
4(6 • 5) 2 (4 • 3) 3 , C 222;22 o = (C 222; ooo) 2 /C2oo;00o — (6 • 5) 2 (4 • 3) 3 , C 222; ooo = 
(6-5)(4-3) 2 , C 22 i ;22 i = 6(6-5) 3 2 5 , C nm2;pg i = C nrn 2- pq o{np+mq + ^) , C nmr . pq o = 

Cnmr;OOoCp<jr;OOo/ Cr00;000; 

v q y 

C 2 22; P q2 = C 2 22; P qi(- + - - 1 + -)• (The coefficient C nmr;0 oo is identical to 



the C nmr of the k=3 case fl29|) . Whenever the graphs reduce to triangles, one 
recovers the coefficients of the k=3 problem). 

However the above representation (|67|) is in fact somewhat ambiguous. 
Indeed there is an interesting cubic identity between the six variables 

Tab = {X a - X b )(X a - X b ) , 1 < a < b < 4, 

which holds for every set of a; a 's and A a 's, namely 

TlV34 + r 2 3 r 2 4 + T 2 4 T 23 + T 23 T U + T 24 T 13 + T^Tu 

— T "l2 r 34(' r 23 + r 24 + r 14 + r 13) — r 13 r 24 ( r 14 + r 12 + r 23 + r 34) 

15 



-^14^23(^12 + T24 + ^34 + T n ) 

+T12T24T14 + 7 "l3 r 14T"34 + T 23 T 24 T 34 + r 12 r 23' r 13 = 0. 



(68) 



Consequently, when one divides the l.h.s. of this identity by the Vander- 
monde squared A(x) 2 A(X) 2 , one generates an identity for the terms of the 
exapnsion of x whose coefficients are C22i ; 220, C22i ; 2ii, Ci22 ; i2i- Thus there is 
an arbitrariness in writing x i n the form of the expansion fl6"7|), in view of 
this cubic identity. For arbitrary k, there are (k — 2)(k — 3)/2 similar cubic 
identities, and therefore only (2k — 3) independent r ab . 

For other values of (3 such as (3 = 8 for instance, by imposing the dif- 
ferential equation (|T9|), one finds that the values of C332 ; 33i, C332 ; 322, C323 ; 322 
are not uniquely determined, because again of this same cubic identity. One 
may need to eliminate this arbitrariness by choosing some arbitrary value for 
some coefficient, decide that is it zero for instance. This cubic identity car- 
ries the problem over to higher orders. Indeed, multiplying this same cubic 
identity, by some symmetric function of the r's, such as VVy, one obtains 
an identity of order r 4 , which again requires some arbitrary choice in the 
expression for x & t next order. 

In summary, we have obtained expressions for the HarishChandra-Itzykson- 
Zuber integration for k x k matrices, for group integrals characterized by the 
parameter (3. Our investigation relies upon new solutions to the heat kernel 
differential equation. Other methods such as supersymmetric quantum me- 
chanics, combinatorial group analysis, or algebraic geometry may be of help 
in solving the numerous open questions concerning this problem. 
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